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We study temperature dependence of geometrical (Fiske) and velocity-matching (Eck) resonances in the flux- 
flow state of small Bi2Sr2CaCu208+x mesa structures. It is shown that the quality factor of resonances is high at 
low T, but rapidly decreases with increasing temperature already at T > 10 K. We also study T-dependencies 
of resonant voltages and the speed of electromagnetic waves (the Swihart velocity). Surprisingly it is observed 
that the Swihart velocity exhibits a flat T-dependence at low T, following T-dependence of the c-axis critical 
current, rather than the expected linear T-dependence of the London penetration depth. Our data indicate that 
self-heating is detrimental for operation of mesas as coherent THz oscillators because it limits the emission 
power via suppression of the quality factor. On the other hand, significant temperature dependence of the 
Swihart velocity allows broad-range tunability of the output frequency. 

PACS numbers: 74.72.Hs, 74.78.Fk, 74.50.+r, 85.25.Cp 
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I. INTRODUCTION 

Single crystals of a cuprate superconductor 
Bi2Sr2CaCu208+x (Bi-2212) represent natural stacks of 
atomic scale intrinsic Josephson junctions (IJJs) LIJ. Joseph- 
son junctions form transmission lines for electromagnetic 
(EM) waves [20. The propagation (Swihart) velocity is 
Co - c/[LaCaV^^, where c is the speed of light in vacuum and 
La and Ca are the inductance and the capacitance per square 
of the transmission line. 



La = 47rA, 
A = t-\-2As coth(d/ As). 



(1) 

(2) 



Here t and d are thicknesses of the dielectric and supercon- 
ducting layers, respectively, and As is the London penetration 
depth of the superconductor. In thin-layer junctions t,d <^ As, 
La - ^ttA^ Id is dominated by a large kinetic inductance of 
superconducting layers. As a consequence, cq can be much 
slower than c - the phenomenon that finds applications in com- 
pact superconducting delay lines [3j . 

The Swihart velocity carries a direct information about the 
London penetration depth. It can be obtained by measur- 
ing the propagation (delay) time in a transmission line [4, 5] . 
However, Josephson junctions provide a much easier way of 
measuring cq. In Josephson junctions EM waves can be gener- 
ated in-situ by means of the ac-Josephson eflTect. At geometri- 
cal resonance conditions they form standing waves, leading to 
appearance of Fiske steps in current- voltage (I-V) character- 
istics f6-l(y\. Fiske step voltages allow siniple and direct eval- 
uation of the absolute values of As (T) i^ luJ] (unlike surface 
impedance measurements, which usually provide only rela- 
tive values [12-15]). Such measurements do not require long 
transmission lines, but can be performed on small ~ yum-scale 
junctions. IJJ of //m sizes, made on high quality Bi-2212 sin- 
gle crystals, are free from crystallographic defects, that can 
aflfect A in cuprates |116i ll?'] . Therefore, Fiske resonances in 
small IJJs should provide information about genuine (defect- 
free) behavior of the penetration depth in cuprates. 



Geometrical resonances play also an important role in 
achieving high power THz EM wave emission from Bi-2212 
mesa structures [118142311 . The maximum radiation power from 
a stack with A/^ junctions is Prad ^ N^Q^ [24], where 



Q = ojRC, 



(3) 



is the quality factor of the resonance, oj is the resonant fre- 
quency, R the eflTective damping resistance and C the capac- 
itance of the junctions. The factor N^ is due to constructive 
interference of N in-phase synchronized junctions [25] and 
the factor Q^ represents the resonant amplification in each 
junction by the geometrical resonance. Thus, both the in- 
phase coherence and the high quality 2 » 1 geometrical res- 
onances are needed for achieving high emission power ll24ll . 
Increment of the emission power is inevitably accompanied 
by self-heating of the stack. In superconductors this leads to 
a rapid increment of the quasiparticle (QP) damping, which 
suppresses Q. Self-heating ultimately limits the performance 
of an oscillator r23] . Clearly, investigation of the quality fac- 
tor of geometrical resonances and their T-dependence has a 
primary significance for development of high power THz os- 
cillator, based on IJJs. 

In this work, we study experimentally T-dependencies of 
geometrical (Fiske) and velocity-matching (Eck) resonances 
[9] in the flux-flow state of small Bi-2212 mesa structures. It 
is observed that Q of resonances is large at low T, but rapidly 
decreases with increasing temperature already at T > 10 K 
<^ Tc ~ 90 K, primarily due to enhancement of the quasipar- 
ticle damping. Surprisingly, it is observed that resonant volt- 
ages, proportional to the Swihart velocity, exhibit a very weak 
T-dependence at low T and do not follow the expected lin- 
ear T-dependence of the eflTective London penetration depth 
Aab(T) in Bi-2212 lll2Ul5l] . We discuss possible origins of 
such a distinct discrepancy, which to our opinion deserves fur- 
ther experimental and theoretical analysis. 



II. GEOMETRICAL RESONANCES IN STACKED 
JOSEPHSON JUNCTIONS 

Stacked Josephson junctions form multilayer transmission 
lines for electromagnetic waves. The general problem of lin- 
ear wave propagation in multilayer transmission lines was first 
considered by Economou |26] and more recently within the 
inductively coupled junction (ICJ) formalism by Kleiner f27'] 
and Sakai et al., 11281] . In this section we will briefly recollect 
peculiarities of wave propagation and geometrical resonances 
in stacked Josephson junctions. 

In the ICJ model of Sakai, Bodin and Pedersen [29], a lay- 
ered superconductor is represented by a stack of isotropic su- 
perconducting layers with the thickness d and the "intrinsic" 
penetration depth As, separated by tunnel barriers with the 
thickness t, the dielectric constant 6^, and the fluctuation-free 
(maximum) Josephson critical current density /cO- The stack- 
ing periodicity s = t+dis^ 1.5 nm for Bi-2212. Properties of 
inductively coupled stacked Josephson junctions are described 
by the coupled sine-Gordon equation [29]. The coupling is 
represented by a tridiagonal coupling matrix A with the of- 
diagonal terms equal to minus the efl'ective inductive coupling 
constant between neighbor junctions [28], 



S =Xs 



t sinh I — 



-h 2As cosh I — 



-1 



(4) 



For atomic scale IJJs, S ^ 0.5 - dslAA^ is very close to its 
maximum value 0.5. 
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FIG. 1. (Color online). Spatial distribution of oscillation ampli- 
tudes of (a) magnetic field and (b) in-plane currents for the in- 
phase (squares) and the out-of-phase (circles) modes for a stack with 
N - \^ IJJs. Horizontal stripes represent superconducting layers. 



is the Josephson plasma frequency and 



Ooc 



87r2/coA 



1/2 



Ooc^ 



len'^JcQX] 



1/2 



(8) 



is the Josephson penetration depth of a single junction and 



A. Eigen-modes in stacked junctions 

The main diff'erence between single and stacked junctions 
is the presence of multiple electromagnetic wave modes in the 
stack. Geometrical resonances in a stack correspond to forma- 
tion of two-dimensional standing waves 123, UM • The wave 
number along the a/7-planes (x-axis) is km = nm/L, where L 
is the length of the junctions and m is the number of nodes in 
the standing wave. In the c-axis direction it is given by one of 
the eigen-modes, kn = 7in/(N -\- l)s,n = 1, 2, ...A/^, where N is 
the number of junctions in the stack. The oscillatory part of 
the phase diff'erence is: 



S(fi(m,n) = a cos 



Inmx 



\ . ( Tini \ 



JCOt 



(5) 



Here / = 1, 2, ...A^ is the junction index, a =const is an ampli- 
tude, and CO is the angular frequency. 

Each eigen-mode has a distinct propagation velocity, given 
by Eq. (3.52) of Ref. ^QM- Within the ICJ model they can be 
written as [28]: 

-1/2 



Cn = Co 



1 - 2S cos 



/ nn Y 

lA^-hi; 



The Swihart velocity cq = Ajojp where 



^TT^tcJr 



O06, 



1/2 



n=l,2,...A^, (6) 



(7) 



^ab - ^s ys/d 



(9) 



is the eff'ective London penetration depth for field perpendic- 
ular to layers. 

Similarly, eigen-modes are characterized by diff'erent char- 
acteristic lengths 130] 



An — A] 



1 - 2^ cos 



/ nm Y 
lA^-hi; 



-1/2 



^=1,2, ...TV. (10) 



(Aj/An)^ are eigenvalues of the coupling matrix A [301]. The 
shortest. An ^ /!// V2 - 0.5 yum for Bi-2212. The longest Ai 
approaches the eff'ective penetration depth for field parallel to 
layers 



Ac = 



Ooc 



87r2X 



cOS 



1/2 



(11) 



for A^ » nAab/s ^ 400. In Bi-2212, AdT = 0) ~ 100 //m 
^ Aab(T = 0)^02 fim^. 

Due to inductive coupling between junctions, the in-plane 
(j-axis) magnetic field is non-local and depends on phase dis- 
tributions in all junctions: By(i) = (Ho/2)A~^Ajd(fij/dx. Here 
Ho = Oo/tt/I/A [30]. Using Eq.© we obtain for the oscilla- 
tory part of magnetic field in the stack: 

HoanmAl 

By(x, z)(m, n) = — sm (kmx) sm (knZ) . (12) 
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FIG. 2. (Color online), (a) Calculated slowest and fastest velocities 
c^ and c\ as a function of the number of junctions in the stack. Cal- 
culations are made within the ICJ model for typical Bi-2212 param- 
eters at r = 0. (b) Calculated temperature dependencies of squares 
of Civ and c\ for several A^ along with /t~^(r), normalized by the cor- 
responding values at r = 0. 



Here we used the property that A and A"^ have the same eigen- 
vectors, and eigenvalues of A"^ are ^J^j, Eq. ([TOl) . 

The in-plane current density in superconducting layers is 
obtained from the Maxwell equation J^ - -{clA7i)dByldz'. 



Jx(x, z)(m, n) = Jacijn, n) sin (kmx) cos (knz) , 

a(!>{)cA^mn 



Jacim, n) 



\6Al^AjL{N+\)' 



(13) 
(14) 



Fig. [T] shows calculated distributions of the amplitudes of 
By (a) and Jx (b) for modes n = N (open circles) and n = I 
(squares) for the stack with N = 10 junctions. Horizontal 
stripes represent superconducting layers. It is seen that the 
eigen-modes are characterized by different symmetry along 
the stacking direction. The slowest n = N mode corresponds 
to the (almost) out-of-phase state in neighbor junctions S(fii ^ 
-S(fi+i . The fastest n = I mode corresponds to the (almost) 
in-phase state S(fi ^ S(fii+i. 

Fig. [2] (a) shows calculated dependence of ci and cn on the 
number of junctions N. It is seen that the slowest velocity is 



almost independent of A^ lliol] 



cn 



Co 
V2 



ts 



^^rAl, 



1/2 



(15) 



To the contrary, the fastest velocity 



ci 




1-F 



^(A^+1) 



-1/2 



(16) 



is growing linearly with N for N < nXabl s - 400 lllOh . 
For A^ » 7iAab(T)/s, it asymptotically approaches the T- 
independent value ci(N -^ oo) = c[t/ErsY^^, close to the 
speed of light in the dielectric, as shown in Fig. [2] (a). 

Fig. [2] (b) shows calculated T-dependencies of c^ and c^, 
Eq. (O, normalized on the corresponding values at T = 0, 
for different N. Calculations are made for typical parameters 
of Bi-2212, using the [^ab(T) / Aab(0)]~^ dependence shown 
by the lowest line obtained from surface impedance measure- 
ments lll2Ul5l ll7'1. As follows from Eq. (O, T-dependence 
of the out-of-phase velocity cn follows l/Aab(T), irrespective 
ofN. For IJJs the same is true for all slow modes n > 2. 

The speed of the fastest mode, ci(T), does depends on N. 
For N < nXabl s - 400, it maintains the same T-dependence 
oc l/Aab' The corresponding three curves [cn(T)/cn(0)]^, 
[Aab(T)/Aabm-^ and [ci(r)/ci(0)]2 for A^ = 100 collapse 
in one in Fig. [2](b). For much larger N, when ci approaches 
T-independent speed of light in the dielectrics, see Fig. [2] (a), 
ci(T) becomes flatter at low T, as shown in Fig. [2](b). How- 
ever, since Aab diverges at T ^ Tc, c\ always vanishes at T^, 
as seen from the curve with N - 10^ in Fig. [21(a). 

In applied in-plane magnetic field Josephson vortices (flux- 
ons) 130,1 enter into the junctions. In strong enough magnetic 
field fluxons form a regular fluxon lattice in a stack. Usu- 
ally a triangular lattice is most stable due to fluxon repulsion. 
However a rectangular lattice can be stabilized via geometri- 
cal confinement is small Bi-2212 mesas 113 ill . Motion of flux- 
ons leads to appearance of the flux-flow (FF) branch in the 
l-V . Emission of EM waves in the FF state leads to excitation 
of geometrical resonances ILZllslllOl]. The corresponding Fiske 
step voltage for the resonant mode (m, ri) is 



ym,niJ) = ^0mCn(T)/2L. 



(17) 



The strongest resonance occurs at the velocity matching 
(VM) condition, when the velocity of fluxons is equal to the 
velocity of electromagnetic waves [lu]. This leads to appear- 
ance of the VM (Eck) step at the end of FF branch ISJ . The 
VM voltage is 



V^ 



VM 



NHsCn. 



(18) 



The T-dependencies of both Fiske and VM steps are deter- 
mined solely by c„(r), Eq.©. Therefore they can be used for 
accurate detection of the absolute value of Aab(T) (except for 
the fastest mode at very large N, as shown in Fig. [2](b)). 



B. Connection between the inductively coupled and the 
Lawrence-Doniach models 



A similar system of coupled sine-Gordon equations was 
also obtained from the Lawrence-Doniach (LD) model 11321] . 
The two main parameters of the LD model are the anisotropy 
factor y =const» 1 and the effective London penetration 
depths Aab- The rest of parameters are derived as [32]: Ac = 



y^ab, ^j = ys, ojp = clsl'^yXab and cq = cslsl'^Xab- 



From comparison with ICJ expressions Eqs. (I7I8I9I1 11151) 
it is seen that while the ICJ model contains two T -dependent 
variables Aab(T) and Jco(T), the LD model has only one, 
Aab(T), which imposes its T-dependence on all other vari- 
ables. Within the range of validity of the LD model, T^ - T ^ 
Tc, the two models are identical because A~^(T) oc Jco(T) oc 
1 - T/Tc. However, as will be discussed below, A~^(T) 
and Jco(T) have distinctly different T-dependencies at low T, 
which does cause a discrepancy between the two models. Es- 
sentially it is related to the fact that in cuprates the anisotropy 
y(T) = Ac(T)/Aab(T) ^const MM- 
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FIG. 3. (Color online) I-V curves of the mesa-1 at // = 1.4 T and 
at different T = 2.0 K - 15.1 K. At low T, in panels (a) and (b), 
sequences of hysteretic (high- 2) individual Fiske steps are seen at 
low bias. At higher bias some junctions switch into the QP state, but 
Fiske steps are still present in the rest of the junctions. The corre- 
sponding first four mixed flux-flow-QP branches are marked (QPl- 
4). At r = 10.1 K (c) these steps smear out and at T = 15.1 K 
(d) individual Fiske steps have vanished, instead a collective, non- 
hysteretic step is observed. 



III. EXPERIMENTAL 

Small mesa structures were fabricated on top of Bi-2212 
single crystals with Tc = 82 K. Twelve mesas with differ- 
ent sizes were fabricated simultaneously on every crystal. 
All of the studied mesas showed similar behavior. Here we 
present data for two mesas on the same slightly underdoped 
Bi-2212 crystal with areas of 2.7 x 1.4 jim^ (mesa-1) and 
2.0 X 1.7 jduprn^ (mesa-2). Both mesas contain N = 12 IJJs. 
The results are representative for a large number of mesas 
made on crystals with different doping and composition (see 
Table-I in Ref.[10]). Details of sample fabrication and of the 
experimental set-up can be found in Ref. lilOl] . 

The magnetic field was applied strictly parallel to the super- 
conducting CuO bilayers, to avoid the intrusion of Abrikosov 
vortices. Eventual entrance of Abrikosov vortices is immedi- 
ately obvious in experiment: it causes very strong and irre- 
versible damping of Fiske resonances and of the Fraunhofer 
modulation of the critical current [31]. Essentially, results re- 
ported here are observable only in the absence of Abrikosov 
vortices. Using the rigorous alignment procedure, described 
in Ref. [31], we were able to p revent Abrikosov vortex en- 
trance in fields up to 17 T ||36i |37|]. This is seen from the 
field-independence of the c-axis QP resistance [36] and per- 
fect reversibility of all measured characteristics [lOl blll. 

All measurements are made in the 3 -probe configuration. 
To simplify data analysis, a contact or a quasiparticle resis- 
tance was subtracted from I-V characteristics, as described in 
Ref. [31]. The subtraction is facilitated by the negligible de- 
pendence of the QP resistance on the in-plane magnetic fields 
due to the extremely large anisotropy of Bi-2212 (see. e.g. 
Fig. 3 (d) in Ref. [36]). To do the subtraction, we first care- 
fully measured the corresponding branch of the I-V at zero 
magnetic field. After that we made a high-order polynomial 
fit of ln(/) vs. y, which is almost linear pS] and can be fitted 
with a very high (~ yuV) accuracy. This fit is then subtracted 
from the measured I-V. When studying T-dependence, this 
procedure was repeated at each T. Such subtraction simplifies 
the analysis of Fiske steps, but is not necessary: Fiske steps 
can be also measured relative to the bias-dependent contact or 
QP voltages. 



IV. RESULTS 

Figure [3] shows I-V curves (digital oscillograms) for the 
mesa-1 Sit H = 1.4 T and at different T from 2.0 K to 15.1 
K. As the current is increased, the I-Vs switch from the zero 
voltage branch to the flux-flow branch, containing sequences 
of individual and collective Fiske steps, seen as small sub- 
branches in Fig. [3] (a), and ending at the velocity-matching 
step. Detailed discussion of the magnetic field dependence 
of Fiske and VM steps at low T can be found in Ref. [10]. 
Strong hysteresis of Fiske steps indicates high g » 1 of the 
geometric resonances. This is facilitated by careful alignment 
of magnetic field, which prevents penetration of Abrikosov 



vortices 13611 . With further increase of current some junctions 
switch into the QP state, while the rest are remaining in the 
flux-flow state. This leads to appearance of combined QP- 
FF families of Fiske steps, four of which are indicated in Fig. 
Ha), (QPl-4) with the number corresponding to the number 
oflJJsintheQPstate. 

The speed can be obtained directly from resonant voltages 
using Eqs.JTTl) and (fTSl) . The corresponding low-T values for 
several mesas at difl'erent Bi-2212 crystals can be found in 
Ref. [10]. Fiske steps in Fig. [3]correspond to slow speed res- 
onances V2,N = 0.27 mV. At the QPl, QP2 branches another 
sequence V4^n = 0.54 mV of individual Fiske steps is seen. 
As shown in Refs. liol [37], the Vvm is proportional to the 
field for 2 T < if < 10 T, consistent with Eq. ([T8]), before 
it gets interrupted by phonon-polariton resonances at higher 
fields 1 37]. In this intermediate field range the limiting fluxon 
velocity is close to the out-of-phase velocity c^. 



A. Temperature dependence of the quality factor 

As seen from Fig. [3l with increasing temperature, the am- 
plitude of the individual Fiske steps rapidly decreases. At 
T = 10.1 K steps are smeared out almost completely and at 
r = 15.1 K they vanish. This indicates a substantial reduction 
of the resonance quality factor. At this temperature only a col- 
lective, non-hysteretic Fiske step is visible at N x V2,n - 3.2 
mV, see Fig [3] (d). 

Figures |4] (a) and (b) show I-Vs in a wider T-range (a) for 
the mesa-2 at ^ = 2.75 T, and (b) for the mesa-1 at if = 3.85 
T. Collective Fiske steps at ^ N x Vi^n can be seen at low T 
(indicated by the downward arrows). At higher bias VM steps 
are observed (indicated by the upward arrows). Both mesas 
show similar behavior: Sharpness of the collective Fiske and 
the VM steps rapidly decreases with increasing temperature. 
This indicates enhancement of damping, also seen from re- 
duction of slopes of I-V curves with increasing T. 

Figure |4] (c) shows dl/dV curves, numerically calculated 
from the I-V curves from (b). Peaks in conductance corre- 
spond to Fiske and VM steps. The decrease of amplitudes of 
the steps with increasing T is clearly seen, indicating reduc- 
tion of Q at higher temperatures. 

According to the sine-Gordon equation, the initial viscous 
part of the flux-flow I-V should be ohmic with the flux-flow 
resistance Rff representing the efl'ective damping 11391] . In- 
deed, from Figs. |4](a) and (b) it is seen that the flux-flow I-V 
is nearly ohmic at 10 < V < 20 mV This allows accurate eval- 
uation of the bare (non-resonant) Rff(T). It is shown in panel 
(d) for y = 12 mV (~ 1 mV per junction). The T-dependence 
of Rff is almost identical to the low-bias c-axis QP resistance 
Rqp(T,H = 0) 11381] . proving that the Rff(T) dependence is 
predominantly determined by "freezing out" of quasiparticles. 
At low T and moderately low H the value of Rff is slightly 
lower than Rqp, which may indicate presence of additional 
damping mechanisms, such as the in-plane QP damping [40], 
or generation of phonons via electros trie tion [37]. At higher 



H, Rff = Rqp (see e.g. Fig. 3 (d) from Ref.f36']). 

Figures in (e) and (f) represent T-dependencies of bare am- 
plitudes of conductance peaks at the collective Fiske step and 
the VM step, respectively. The peak amplitudes were ob- 
tained by subtracting the background flux-flow conductance 
R~j}p. It is seen that resonances in both mesas exhibit similar T- 
dependencies: At low T, peaks are high, i.e., quality factors of 
resonances are large 2 » 1 , but they start to rapidly decrease 
with increasing T . Comparison with the efl'ective flux-flow 
resistance Rff, shown in panel (d), indicates that the scale for 
variation of peak amplitudes is similar to Rff(T). Therefore, 
both resonances roughly follow Eq. Q with R ^ Rff(T). 



B. Temperature dependence of the Swihart velocity 

Both Fiske and VM steps in the considered case correspond 
to propagation of waves, respectively fluxons, with the veloc- 
ity ^ 3.2 X 10^ m/s 1 10] is close to the expected value of the 
slowest out-of-phase velocity cn, Eq. ([T5]) . It is almost 1000 
times slower than c, not because of extraordinary large di- 
electric constant, but because of extraordinary large kinetic 
inductance of atomically thin superconducting layers in Bi- 
2212, see Eq. (2). According to Eq. ([T5]) , cn(T) should 
depend solely on l/Aab(T). Thus voltages of Fiske and VM 
steps should provide a direct information on absolute values 
of I /Aabin 

Squares and triangles in Fig. [5] represent measured T- 
dependencies of Vyj^ for both studied mesas. Crosses in Fig. 
[5] represent fast geometrical resonance voltages, reported re- 
cently by Benseman and co-workers on large Bi-2212 mesas 
at zero field |41]. Apparently, our data for the slowest reso- 
nances coincide with their data for the fast resonance within 
the measured T-range. 

Lines in Fig. [5] represent typical temperature dependencies 
ofA~^ for cuprates Ill2l[l3|] and the fluctuation-free c-axis crit- 
ical current density /^o for Bi-2212 IJJs |1, 42]. The latter 
is similar to 0J^(T), measured by the Josephson plasma reso- 
nance 11431] and to A~^(T) obtained from surface impedance 
measurements ll34|] . consistent with Eqs. (l7llll) . It is seen 
that AZu and Jco exhibit distinctly difl'erent behavior at low T: 



ab 



Jco(T) is flat, while A~^(T) has a linear T-dependence due to 
the d-wave symmetry of the order parameter Ill2lll3|]. Clearly, 
experimental Vyj^(T) follow Jco(T) rather than the expected 
^ab(T)~^ dependence. 



V. DISCUSSION 

At low r, the obtained speed of EM waves ^ 3.2 x 10^ 
m/s agrees with the expected out-of-phase mode velocity cn, 
Eq. ([T5]) for reasonable para meters t/Sr = 0.1 nm and 
Aab(T = 0) ^ 200 nm O El Hill. Thus the ICJ model 
does provide a correct value of the Swihart velocity at low 
T. It also provides correct T-dependencies of the Joseph- 
son plasma frequency [.43.] and Ac 134], ojp(T) oc A~^(T) oc 
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FIG. 4. (Color online) (a) Flux-flow parts of I-V curves of mesa-2 at // = 2.75 T and at diff'erent T. A collective Fiske step (downward arrow) 
and a velocity matching step (upward arrow) are seen, followed by QP branches at higher bias V > 20 mV. (b) The same for the mesa-1 at 
// = 3.85 T. It is seen that both Fiske and velocity matching steps are rapidly smeared out with increasing T. (c) dl/dV curves, numerically 
calculated from curves in (b). Distinct peaks correspond to the collective Fiske and the velocity-matching resonances, (d) T- dependence of the 
nearly ohmic flux-flow resistance at V = 12 mV for the mesa-1 at // = 3.85 T. (e) and (f) show amplitudes of the dl/dV peaks, corresponding 
to the collective Fiske step (e) and the velocity-matching peak (f). 



V/co(^), see Eqs. (l7llll) . Therefore, it is surprising that the T- 
dependence of the effective penetration depth deduced from 
resonant voltages is different from Aah(T), obtained from sur- 



face impedance measurements |112l [Ij, [151 133|]. Below we 






mention several possible reasons for such a discrepancy. 



A. Possible origin of discrepancy with surface impedance 
measurements 



Derivation of the ICJ model is based on the assumption 
that field and current distributions within each superconduct- 
ing layer can be described by the local 2nd London equation 
II29II . However, this assumption most likely breaks down in 
atomic scale IJJs (see the condition (4.3) in Ref. |26]). 

To understand the reported discrepancy it is, first of all, nec- 
essary to understand the difference in local current and field 
distributions. In surface impedance measurements, the exter- 
nal electromagnetic field is screened at the depth Aab ~ 200 
nm from the surface of the superconductor. This induces 
similar (in-phase) screening currents in a fairly large number 
N - 130 of IJJs. To the contrary, at the out-of-phase geo- 



metrical resonances the current varies at the atomic scale, as 
shown in Fig. [l](b). 

/. Non-locality of supercurrent 

The most obvious question is to what extent Cooper pairs 
are localized in every CuO bi-layer. The very existence of the 
c-axis critical current indicates that the localization is incom- 
plete. This can be particularly significant for the out-of-phase 
mode, when Cooper pairs are forced to move in opposite di- 
rections in neighbor layers, see Fig. [l](b). Qualitatively such 
delocalization will lead to larger effective penetration depth. 

//. Nonlocal Josephson electrodynamics 

Another type of non-locality in thin layer junctions was 
considered in Ref. IA4\. With decreasing d, the effective 
screening length A/ 2 Eq.(3) increases and approaches the 
Pearl length Ap = A^/d. To the contrary, the Josephson 
penetration depth Aj decreases oc hr^l^ see Eq.®. For IJJs 
Aj < I jum [31] is much smaller than Ap - 10 yum even at 
r = 0. Such a mismatch changes the dispersion relation of 
electromagnetic waves 1441 . 

///. Retardation effects 

Retardation effects appear in transmission lines when the 
time (phase velocity) required to transfer charge within a layer 
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FIG. 5. (Color online) Normalized temperature dependence of the 
square of the velocity-matching voltage Vyu ^ c^- Solid line repre- 
sents typical T-dependence of A~^ from Ref. |12, 13]. Dashed line 
represents the T-dependence of the fluctuation free Josephson criti- 
cal current density 7^0 ^ ^l from Ref. 11421] . 



is comparable or faster than that for electromagnetic waves 
outside the layer [,261 . Specific for IJJs is that the out-of-phase 
electromagnetic wave velocity is so slow ~ 10^ m/s, that it 
becomes comparable to the electronic Fermi velocity. This 
may aff'ect the dispersion relation. 

iv. Frequency dependence 

The eff'ective penetration depth in superconductors depends 
not only on T but also on frequency A{T, oj). It originates from 
a significant (T, o^ dependence of complex conductivity in a 
superconductor |11]. The most obvious diff'erence between 
static and high-frequency A(T) is that the latter does not di- 
verge SitT ^ Tc, but approaches the finite normal skin-depth. 
This may flatten-out T-dependence of high frequency Fiske 
steps, compared to static l/A(T) |11]. 

Surface impedance measurements are typically performed 
at ~ 10 GHz frequency. In comparison, the studied Fiske 
and VM step voltages are ~ 1 mV per junction, see Fig. |4l 
According to the ac-Josephson relation, this corresponds to 
-500 GHz. The significant diff'erence in frequencies may 
lead to a significant diff'erence in the eff'ective A. 

At even higher THz frequencies, the frequency dependence 
of the dielectric function Sr(oj) in isolating BiO layers be- 
comes significant. As shown in Ref. |37], the speed of electro- 
magnetic waves slows down dramatically, when the frequency 
approaches the transverse optical phonon frequencies. 

V. Non- linear effects 

Eq.® was derived by linearization of the coupled sine- 
Gordon equation and is valid for small amplitude EM waves 
a <^ I. However, at high quality geometrical resonances the 
amplitude may be large a - 27i and non-linearity of the sine- 
Gordon equation may aff'ect the dispersion relation. 

The penetration depth depends on the absolute value of 



the current density. Close to the depairing current density, A 
rapidly increases. At geometrical resonances the amplitude of 
the in-plane current density is given by Eq. (14). It depends 
on the amplitude a, which can be ~ 1 for 2 » 1. An esti- 
mation for a = 1 and L/m = 1 yum yields Jac -10^ A/cm^, 
comparable to the maximum in-plane current density 11451] . 

Both types of non-linear eff'ects increase with increasing the 
quality factor of resonances. Since 2 » 1 only at low T, 
nonlinear corrections can be significant at low T, but less so 
at elevated temperatures. 



B. Implications for coherent Josephson oscillators 

As mentioned in the Introduction, stacked IJJs are consid- 
ered as possible candidates for high power THz oscillators 
Ili-IIlEil]. A large energy gap in Bi-2212 fSIH allows 
generation of electromagnetic radiation with frequencies in 
excess of 10 THz. For example, recently polariton gener- 
ation with frequencies up to ~ 13 THz was reported |37]. 
Moreover, strong electromagnetic coupling of IJJs facilitates 
phase-locking of many junctions, which may lead to coherent 
amplification of the emission power 1,25,1 . 

Realization of a flux-flow oscillator yjL based on fluxon 
motion in the in-plane magnetic field [^ I46l448l] . encounter a 



difficulty, associated with instability of the rectangular fluxon 
lattice. It can be stabilized by geometrical confinement 
in small mesas 113 ill or by interaction with infrared optical 
phonons pT"]. But usually fluxon-fluxon repulsion promotes 
the triangular fluxon lattice, corresponding to the out-of-phase 
state, which leads to destructive interference and negligible 
emission [24]. 

High quality geometrical resonances improve operation of 
a stacked oscillator is several ways: (i) they amplify the emis- 
sion power oc Q^ [24]; (ii) they narrow the radiation linewidth 
oc 1/2 [24]; (iii) they C3.n force phase-locking of junctions. 
Numerical simulations have demonstrated that large ampli- 
tude standing waves, a ~ 1, can superimpose their symmetry 
on the fluxon lattice 1 10]. Such a non-linear synchronization 
requires high Q because a oc l/Q. 

The reported rapid decrease of the quality factor with in- 
creasing temperature indicates that self-heating is detrimental 
for the coherent Josephson oscillator and ultimately limits the 
emission power from large Bi-2212 mesas |23]. On the other 
hand, T-dependence of the Swihart velocity facilitates fairly 
broad-range tuning of the resonance frequency, as seen from 
Fig. \5\ This may be beneficial for the oscillator 1.4 1.1 . 



CONCLUSIONS 

To conclude, we have studied T-dependence of geometri- 
cal and velocity matching resonances in small Bi-2212 mesa 
structures. We reported strong T-dependence of the qual- 
ity factors, which is large at low T, but rapidly decreases 
with increasing T, already at T > 10 K. Above T ~ 60 K 



~ 0.87^ resonances are almost fully damped. This observa- 
tion is consistent with previous observations of strongly un- 
derdamped phase dynamics at low T 143], leading to rela- 
tively high macroscopic quantum tunnelling temperature in 
IJJs 1I49I451I1 . and with the reported collapse into overdamped 
dynamics at T jTc -0.8 11421] . The rapid decrease of Q{T) 
indicates, that self-heating is detrimental for operation of the 
coherent THz oscillator and ultimately limits its performance 
1I23I] . On the other hand, T-dependence of the S wihart veloc- 
ity facilitates a broad-range tuning of the resonance frequency, 
which may be beneficial for the oscillator. 

Our analysis of T-dependence of resonant voltages revealed 
that the eff'ective penetration depth, that determines the ki- 
netic inductance and the speed of electromagnetic waves in 
intrinsic Josephson junctions, see Eq. ([B, is exhibiting a flat 
T-dependence at low, resembling T-dependence of the c-axis 
critical current. It is distinctly difl'erent from the linear T- 
dependence of A^(T), obtained from surface impedance mea- 
surements III2I ll3[] . We argued that non-trivial physical phe- 
nomena, such as break-down of the local London approxima- 
tion at the atomic scale, are responsible for this distinct dis- 
crepancy, which deserves further theoretical consideration. 
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